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Abstract 

The TexSol is a composite geomaterial : a sand matrix and a wire network re- 
inforcement. For small strains a thermodynamical continuous model of the TexSol 
including the unilaterality of the wire network is postulated. This model is described 
by two potentials which depend on some internal variables and a state variable either 
strain or stress tensor (the choice of this last one gives two different ways of identi- 
fication). The TexSol continuous model is implemented in a finite element code to 
recover the mechanical behaviour given by discrete elements numerical experiments. 
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1 Motivations 



1.1 What IS the TexSol ? 



The civil pieces of work need planed stable floor. The environment configura- 
tion often forces civil engineers to raise huge embankments. Moreover, it can 
be interesting to reinforce them to assure a better embankment mechanical 
behaviour. A lot of different solutions can be used to reinforce soil but, in this 
paper, we focus our attention to the TexSol process. 

The TexSol, created in 1984 by Lefiaive Khay and Blivet from LCPC (Labo- 
ratoire Central des Fonts et Chaussees) [12], is a heterogeneous material by 
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mixing sand and wire network. This reinforced material has a better mechan- 
ical resistance than the sand without wire. Of course, the TexSol behaviour 
depends on sand and wire parameters and its frictional angle can be larger 
than the sand one from 0° to 10° [8]. The wire is described by its linear den- 
sity with a dtex unit (1 dtex = 0,1 g.km~^), its ponderal content and its 
stiffness. Classically, the wire density in a TexSol sample is included between 
100 km.m~^ and 200 km.m~^. 

To make a TexSol bank, a machine named "Texsoleuse" is used. It works on 
throwing sand and, in the same time, injecting wire. The wire is deposed on 
the free plane of the sand with a random orientation. This machine carries 
out several passes to raise the bank. The figure 1 is the TexSol microstructure 
representation. In the literature, we find two different continuous modellings. 




Fig. 1. Schematic TexSol sections 



The model suggested in [5] is non local and includes remote interactions (cor- 
responding to the wire effects) but also needs an identification of their param- 
eters using macroscopic experiments. Villard proposes a simpler local model 
in [20]. This one couples a standard model of sand and an equivalent unilat- 
eral elastic stiffness contribution corresponding to the wire network. This last 
contribution is activated only on the traction directions because of the unilat- 
eral behaviour of wire. Our first work is to clearly define thermo-dynamical 
potentials of the Villard local model with both stress and strain formulations 
to identify the best-adapted one. 



1.2 Assumptions of the continuous local model 

To couple the elastic plastic model of the sand and the unilateral elastic model 
of the wire network, we have to consider some mechanical assumptions, which 
may be backed up by numerical experiments performed with a discrete ele- 
ments software [3,15]. 
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1.2.1 Stress additivity assumption 

In this paper, the stress additivity assumption of the sand and the wire network 
is assumed. Then we write, 

(Ts + (Ty, = a- , (1) 



where cr^, cr^ and cr are respectively the stress second order tensor of the 
sand, the wire network and the TexSol. 

This assumption seems to be coherent with the TexSol quasi-static behaviour. 
We can get a good approximation of the stress tensor in numerical simulation 
of 2D granular matter [16] using the Weber stress tensor [1]. This tensor 
may be non symmetrical if inertial effects are not negligeable. For quasi-static 
processes this discrete tensor is a good candidate to represent a continuous 
stress tensor. Moreover we can define such a tensor grain by grain with the 
Moreau approach [15]. In this way a wire network stress and a sand stress may 
be computed, to recover by addition the full TexSol stress (in the simulation, 
the wire is modeled by a chain of beads with unilateral interactions [11]). 
On a biaxial crushing test we verify the symmetry property even for large 
deformation as long as the process remains slow. 




-l,2e+05 



Fig. 2. Stress eigen values evolution in the TexSol 
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The eigenvalues (EVl, EV2) are computed and the contribution of each com- 
ponent of the TexSol are underhned in the figure 2 where the wire network is 
only in a tensile state ; in the two eigen directions the sand is in compression. 
This may be also observed on the distribution of force network in the granular 
sample. 

1.2.2 Non sliding assumption 

This second assumption is not as evident as the previous one. Although micro 
slidings occur between sand grains and wire, we assume that at the macro- 
scopic level of the continuum model, the sand network does not slip through 
the wire network. This assumption can be translated by the equahty of the 
three strain rates, 

£5 = = e , (2) 

where e^, £^ and £ are respectively the strain second order tensor of the sand, 
the wire network and the TexSol. 

We have to be very careful with such a condition and define some validity 
domains for it. Indeed the limits of this assumption are difficult to quantify 
and we will restrict the validation of the following continuum model to small 
strains. 

1.3 Role of the wire unilaterality 

The wire network contributes to the tensile srtength of the composite mate- 
rial but not to the compression one (cf. figure 2). To model such a behaviour 
at the macroscopic scale, it is convenient to introduce a unilateral condition 
in the behaviour law of the wire network. This unilaterality accounts for two 
microscopic phenomena. The first one is the lack of bending strength of the 
wire network viewed as a piece of cotton. The second one is the local buckling 
of short segments. The first aspect is not explicitely taken into account by a 
unilateral condition at the microscopic scale in our discrete numerical simula- 
tion since the chain of beads has no bending strength. The second aspect may 
be enforced by introducing a unilateral interaction law between two successive 
beads. Such an interaction models an elemental wire between two beads : we 
denote by "rigid wire"; if not we speak about "rigid rod" (cf. figure 4) for 
bilateral law between beads. The figure 3 illustrates the difference of global 
behaviour between both simulation for crushing biaxial test. Until 6 percents 
of deformation, the responses are almost identical. For larger deformation the 
"rigid rod" model leads to a rough increase of rigidity due to the appearance of 
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Fig. 3. Force / Crushing response with different models 
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Fig. 4. Schematic discrete laws 

compression columns in the wire. Such a phenomenon seems not very reahstic 
and is probably issued from a scale effect since the numerical sample is not 
representative enough of the material. In particular the model of wire with a 
chain of beads generates non realistic wedges of beads by sand grains. 



l.J^ Why a strain and stress formulation ? 



In this paper, we propose to carry out a thermodynamical study with both 
strain and stress formulations. The interest of this work is in the identification 
possibilities of potentials parameters. Indeed, an experimentator making some 
tests on a sample has only access to the global strain. Our numerical inves- 
tigations allow us to have access to finer data such as the local stress field 
throughout the sample. Moreover the global stress tensor over the sample can 
be deduced by an average. 

The post processing of numerical experiments mentioned in §1.2.1 provides 
precise informations on the stress fields, in the sand and in the wire network. 
The stress "unilaterality" in the wire is clearly established in the figure 2. This 
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observation could lead us to favour a stress formulation. But the finite element 
softwares are essentially developed using a strain formulation. Consequently 
we propose, in the following study, strain formulations easily implcmentable. 
Dual stress formulations are provided when they can be analytically deduced 
by the Legendre - Fenchel transformation. 



2 A general thermodynamical framework 

In this part, we define potentials written with different state variables. These 
potentials have to check the Clausius-Duhem inequality to be thermodynam- 
ically admissible. 



2.1 Strain versus stress approach in thermodynamics 

This work must be as exhaustive as possible, while passing from unspecified 
state variables to its dual. We thus use the Lcgcndre-Fenchel transformation 
[14], to carry out our study with both strain and stress formulations. Let us 
write the Clausius-Duhem inequality where u is the internal energy, s the 
entropy, q the heat fiow and T the temperature. 



The intrinsic dissipation depends on a state variable X (or its dual X*), 
some internal variables a = {ai, ■ • ■ , ap} (each internal variable can be scalar, 
vectorial or tensorial) and the temperature T. It can also be expressed with 
the free energy t/j or its Legendre-Fenchel transformation ip* with respect to 
the state variable X, 



where X is the argument of the supremum. Considering either X or X*, we 
find two expressions of the Clausius-Duhem inequality. 




(3) 



i;{X, a, T) = u{X, a, T) - Ts{X, a, T) 
and r{X\a,T)^s^{X:X*-i;{X,a,T\ 



X 



')} 



= X : X* -il){X,a,T) , 



a -.£ — 



(4) 
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(5) 



Using the Helmholtz postulate (this last one can be applied with the gener- 
alized standard materials assumption [6]) and the previous definitions, we are 
now able to deduce the state laws, 



Primal state laws 



Dual state laws 



X* e dxi'iX,a,T) 
—s G dxi^^X, a, T) 



X e ^x*^P*{X*,a,T) 
s e dTiJ*{X*,a,T) 



(6) 



where is the thcrmodynamical force associated with a^. Formally we use 
subdifi^erentials instead of derivatives. If convexity is not required, previous 
relations still hold using the Clarke subdifferential [2]. Then the primal and 
dual forms are not necessary equivalent. In the general case, the Clausius- 
Duhem inequahty (4) or (5) can be reduced to a dot product of a vector flow 
and a vector force. 
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where 



e or e 
cr*'' or cr 



(7) 



The flow variables have to be related by evolution laws to the force variables. 
To verify the inequality (7) some assumptions may be added to these relations. 

It is convenient to introduce a dissipation potential (/? from which the evolution 
laws are derived. By duality a force function (p* is automatically defined using 
the Legendre-Fenchel transformation. 



Primal complementary laws 



Dual complementary laws 



5- G d^ip{e, a, 

£ 

Aed^^{e,a,-^) 



(8) 
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To verify the Clausius-Duhem inequality, some assumptions on the dissipation 
potential are necessary. For simplicity we consider now an isothermal process. 
The left-hand side of the inequality is reduced to, 



and the primal state laws are summarized in (ct, A) G dLp{£,a). ip being a 
separately convex function, with a convex analysis characterisation of the sub- 
differential we write. 



Moreover, if ip is minimum in (0,0), the Clausius-Duhem inequality is then 



Similar properties are required for ip* to recover the Clausius-Duhem inequal- 
ity. Generally we distinguish the reversible and irreversible parts of the trans- 
formation. We thus postulate an additive decomposition for both reversible 

and irreversible parts of the strain tensor e = + e**" and the stress tensor 
cr = cr^' + cr'''^. The reversible / irreversible splitting of cr is less classical. To 
illustrate its interest, remark that eventual residual stresses may be accounted 
for in the irreversible part. 

At this stage we have to choose the external state variable X for the strain 
formulation and consequently X* for the stress formulation. It is usual to 
consider for X the total strain tensor e. By the way the reversible stress <t^ 
appears in the state law and becomes the state variable in the dual stress 
formulation. But we can also use the reversible strain part and deduce the 
full stress tensor cr as the dual state variable (cf. table 1). The first column 
expresses the primal model using e** or e as state variable. The second one 
provides the corresponding dual formulations. 

2.2 ID model of reinforced geomaterial 

Let us apply previous results to a rheological ID model of TexSol taking into 
account the wire unilaterality. 

2.2.1 Strain formulation 

We choose to superpose a classical ID model of elasto-plasticity with hardening 
for sand [13] and a ID unilateral model of elasticity for wire. We thus propose 





verified [18], 



a -.e+Aay Lp{£, a) - p{0, 0) > . 
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Table 1 

Strain versus stress formulations 

the two potentials (free energy) and (/? (dissipation potential) depending 
on the external state variable e and on the internal one £2 as sho-wn in the 
figure 5, 
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Fig. 5. Rheological TexSol diagram 



'il)i{e,e2) if £ e Ci 
■02 (£,£2) if £ e C2 



(9) 



ip{e, £2) = ctq 



£2 



(10) 
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where , and 

^2(£,£2) = \Ks{e-e2f + \He2^ C2 = G R|£ < 0} 

(7o the stress threshold. According to the table 1 the state and complementary 
laws are derived, 

State laws 

(7^ = f = K^{e) + K,{e-e2) , A = -g = {e - £2) - He2 

Complementary laws 
(7^^ = ^ = , A e 5 . 

de £2 

where (e) = max(0, e), the non negative part of e. 



\^sign{e2)(To^ if £2^ R* 
[-(To, (7o] if £2= 



^.^.^ Stress formulation 

To determine the stress formulation, we have to calculate the Legendre- 
Fenchel transformations of ip and (p which are not always analytically accessi- 
ble. However we can use the following general result convenient for piecewise 
smooth functions. 

Proposition 1 Consider a non overlapping splitting (Ci),-^^ ^ of the strain 
space R'^^'^, U7=i = R"^^'^, Q close convex cone with mes{Ci fl Cj) = 0, 
i 7^ j. Ift(j{£) is piecewise defined by ■0(e) = if £ £ Ci, i — 1, . . . ,n then 



ij*{a) = sup{[riV'^c°) (cr)} 



Proof : Let recall the definition of inf-convolution of two functions / and g 
[14], the indicator function of a convex set A and the polar cone C° of C, 



o (/V^)(^)= inf {/(£i)+/(£2)} 

e=ei+e2 



if £ e A 
+00 if £ ^ A 
o C° = {<r |e : <r < , Ve e C} 
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According to classical rules of convex analysis, 



ip*{a) = supjfT : £-inf{V'i(e) +Ici(£)}| 
= sup jsup {<T : e - ^^(e) - Ici(£)} j 
^ sup {{^|;, + lc,n<T)} 

i 

= sup{(^* vlc°) (o-)} . 



For the ID model the splitting into two half spaces is obvious and the analytical 
forms of conjugate functions from (9) are reachable, 



2Ks 2 



, C° = {(T e R|(j > 0} . 



Using the proposition 1, we obtain successively. 



-02 V Ic° = < 



if o-'-+A's£2<0 



2{K^+Ks) 2 " +K,e2>0 



(Ks+H)s2^ 



if (7-^+Kse2>0 



(a^+K,e2f _ (Ks+H)e2^ -r 

2 +Ks£2<0 



Finally, 



^ (^:±M! - i^^i^l^ if + < 



fill 



The Legendre-Fenchel transformation of the dissipation potential is computed 
classically from (10), 



^*K,£2)=I{0}(0-^0 



£2 



(12) 



We implicitely get from the equation (12) : a*'' = 0. The state and comple- 
mentary laws in the stress formulation are straightforward derived. 
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State laws 



if a'-+KsS2>0 



if <t''+Xs£2<0 
K^VV, [T''+Kse2)-[Ks+H)e2 if (t'-+Js:.£2>0 



a'-+Kse2-{Ks+H)e2 



iia'-+Kse2<0 



Complementary laws 
_a.^*((T''-,£2) 



{|sifirn(£2)o"o} si R* 
[-(To, do] si £2= 



This set of equations is equivalent to the one obtained with the strain formu- 
lation §2.2.1. 



3 Strain and stress approach for 3D models 



The complex microstructure of the TexSol material needs not to neglect the 
three dimensional effects. To define a 3D model we follow the previous ID 

approach superposing a classical elastic plastic behaviour for the sand and a 
unilateral elastic one for the wire network. Simple and sophisticated unilater- 
ality conditions may be considered leading to different formulations more or 
less easy to handle in a general primal / dual framework. 



3. 1 3D thermodynamical potentials of the sand 



First of all, let us recall that the stress tensor can be split into a spherical part 
and a deviatoric one, 

<T = Jtr(<r) J + S 

"——^ ' deviatoric 

spherical 

Let us introduce the spherical projection tensor S = and the deviatoric 

projection tensor ID = I — S. In a classical model the state variable is the sand 
full strain £g, the internal ones contain the plastic strain e^, the kinematic and 
isotropic hardening variables a and p [21]. The free energy ipg has the following 
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form, 



£f , a,p) = ^ (e, - : K, (s, - e^) + : a + ^^9^ . (13) 



where K^, Hk and ff^ are stiffness coefficients. Tfie state laws are directly 
derived from it. 





des 




A = 






X = 


dills 
da 




R = 


dtps 
dp 





(14) 



To derive the complementary laws it is more convenient to define the force 
function (p* instead of the dissipation potential (p, 

(^:«, A, X, R) = I{o}«) + In(x,«)(A) , (15) 



where Q{X, R) = {A \F{A, X, i?) < 0} is the elastic domain bounded by the 
Drucker - Prager criterion F defined by [4], 



F{A, X, R) = ^Jj2{A - X) - ry{A) - R{p) . (16) 



Remark that y J2{-) is the pscudo norm of the tensor deviatoric part implied 
in the plastic phenomenon. The initial threshold Ty depends on the pressure 
(as it is usual in soil mechanics), on the friction coefficient P related to the 
friction angle and on the cohesion parameter Co, Ty{A) = Co — /3tr(A) = 
Since we use the dual dissipation potential, we get the complementary laws 
usally issued from the stress formulation (cf. table 1), 



eed^ir^*{(Ti',A,X,R)^B:' 



3x3 



P =X%(A,X,R) = -A , 



A-X 



(17) 
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where A is the plastic multipher always non negative. Its value can be found 
with the plastic condition F — and the consistance condition F— 0, 



F = 
F= 



MA-X)^Ty{A) + R{p) ^^^^ 

, (2(r,(l)+fl(p)) + (^^) ■ 



Contrary to the ID case, we cannot explicitely express a 3D dissipation po- 
tential depending on flow variables. 



3.2 Unilateral wire network model 



According to the ID model, we neglect the dissipation effects, and we focus on 
the free energy. Its stiffness cannot be reduced to the stiffness of the wire and 
has to account for the wire distribution in the sample, assumed to be isotropic 
in the following. Due to the entanglement of the wire network, it is convenient 
to consider continuously differentiable free energy to derive smooth relations 
between strain and stress at the macroscopic level. A model directly derived 
from the isotropic linear elasticity may be expressed in the eigen directions ; 
the strain and stress tensors have the same ones. Consequently the free energy 
is simply written using the Lame coefficients X^, and the strain eigen 
values denoted e^, (we introduce the notations = diag(£:^, e^, e^) 

and (e^) = diag((4), {el), (4))), 

Me.) = + 4 + + ((4)2 + (4)2 + (4)2) . (19) 



The first term describes the volumic unilateral behaviour of the wire network 
activated by the trace of the strain. The second part concerns the shear com- 
ponent which is not activated in all directions simultaneously but according to 
the sign of the strain eigen values. The stress expressed in the eigen directions 
is easily derived from this previous energy, 

a-l = ^t^isw) = A^(tr(e^))I + 2//^(e^) . 

In the current frame the strain stress relationship has the form, 

< = A^(tr(e^))/ + 2ii^P{e^)P^ , (20) 

where P depending on is the passing matrix from the eigen directions to 
the current ones. The expression P(i^)P^ is called the positive part of the 
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wire strain tensor denoted e^. The convexity of the free energy is an open 
question in the general case but it is easily verified for jj,^ — because the 
trace is a linear operator. 



3.3 Models superposition and TexSol potentials 



The previous model is combined according to the ID approach. Moreover, 
we introduce two eventual initial stresses cr° and cr°. There are generated 
by the deposit process under gravity which may be simulated by a discrete 
element software [3,15]. Then we can reasonably assume that eigen values of 
<T^ are non negatives. We define the corresponding the initial strains using the 
elastic parts of the previous models, £^ = K~^cr° and £° = Kj^cr°, where 
— ^vjl ® I -\- 2/Xt„I. The total free energy is then postulated, 

^(£,£f,a,p) = ^(e-e^' + eO) :K,(e-eP + £2) 

+ ^(tr(s + el)f + (s + slf : (s + slf (21) 
: a + ^p^ . 



The state laws are derived, 

o-^ = f = K, (s - en + cT^s + A«,(tr(£ + el))I + 2//^ (e + elf 



The complementary laws are derived considering the dual dissipation poten- 
tial of the sand alone (cf. equation (15)). In the simple case where yu^ = 0, 
we can complete the dual stress formulation by computing the Legendre - 
Fenchel transformation of the free energy ip (denoted in this case 'i/'o) via the 
proposition 1. 



i(o-'-+K,(ef-£2)+o-0):(K,+K°)-i(o-'-+K,(£f-e2)+<TO) 

-i((ef-£«):Ka(eP-eO)+/ffca:a+J/ip2-)_i^o^^o 



V'o(o•^£^,a,^>) 



if tr(cr^ + K, {eP-e1-el)) > 

i (o-'-+K, (ef -eO)) :K,-i (o-'-+K, (e^-e^)) 
-i((£f-e2):K.(ef-eO)+iffca:a+i?ip2-) 

iftr(o-'- + K, (ef-£0-£0)) <0, 



(23) 
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where K° = A^J I. The Legendre - Fenchel transformation cannot be 
catched in the more general case. 



4 Numerical development 

Starting from a coherent thermodynamical model for the TexSol, the next step 
consists in implementing it in a finite element software [7,9]. We discuss then 
responses provided by the simulation of simple compression / traction tests 
according to the expected behaviours detailed in section 1. 



4-1 Numerical implementation 

The variables being known at step n — 1, we have to compute them at step 
n using a predicted value of the strain increment A£„. In a sake of simplic- 
ity, the initial stresses are neglected (cr^ = cr° = 0). Two sets of variables, 
(cr5 „,X„,p„) for the sand and (cr^ ^) for the wire network, are computed si- 
multaneously. The stress in the wire network is directly deduced from the 
potential defined by (19). For the sand the relations given in (14), (17) and 
(18) can be reduced to three equations depending on the three unknowns 
(cTg „, X„,p„). This system is solved by a Newton - Raphson method applied 
to the following residuals ; a = 1,2,3. 

— ^ (cr s,n — Cr s,n-l) 

Q2 Hi;{p„—pn-i) ( c v^-LV y 

Qn — {Ra+Ty^„){Hk+2Hi) i^s,n " + 2 (i?„ -|- %n) PI) '■ {<^s,n ~ (^s,n-l) 
+Pn -Pn-1 , 

where a — 1 corresponds to equations (2), (14)i, (17)2,4 and (18)i, a — 2 
corresponds to equations (14)3, (17)3,4 and (18)i and finally a = 3 corresponds 
to equations (17)4 and (18)i,2 (in all these equations, i?„ is calculated using 
the equation (14)4). Classically, the Taylor development is defined as follow. 

The analytical formulations of the tangent matrix coefficients are given. 
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[{Rn + ry,n) ID + /3 {Ss,n - Xn) /] - K, 



dQi _ Pn-Pn-1 TT 

_ -Rn+ra,n+-gi(p»-Pn-l) ( Q "V \ J_ T 

l2 



(1 - Hk^^) I 

TT Rn+Ty,ri+Hi{p„-pn-l) / C" _ Y 



9Q: 



'^^n {^s,n ~ '^n) • (f s,n ~ s,n-l) I 
QX ^71 \ s,n '-'s,n—lj 

^ = 1 + Cl {Ss,n -Xn-2{R^ + Ty,n) 131) : {,<7s,n ' <y s,n-l) , 



where Ci(^n,rj,,„) = (r„+^^,„)V,+2H.) ' ^n(^n,rj,,„) = ^^f^ and finally 

{Rn,Ty,n) = j^'^^^" • The algorithm Is schematized In the table ?? (where 

(Cn) Cn; Cn) = (<''s,n, ^n, Pn)) ■ Thls last One being quite complex for the sand, we 
have compared the results given by the previous Integration law and strategy 
with the one developped In the CastSM software where a Drucker - Prager 
finite clement model Is available. Since we got a good agreement with both 
Implementations, we focus our attention on the coupled sand/wire model of 
TexSol Involving a unilateral behaviour. 



4.2 Patch test 



In a first step, the simple patch test considered Is a single Ql-Lagrange hexahe- 
dron finite clement submitted to a traction/compression loading (cf. figure 6). 
More precisely, a confinement pressure Is prescribed via a cohesion behaviour 
on the material [17] depending on a single coefficient Cq. A displacement Is 
Imposed on the upperslde. Four models are compared to underline the per- 
tinency of the two unilateral behaviour laws. Two of them are considered to 
obtain some limit behaviours ; the first one denoted Sand, Is free of wire ; 
the second one denoted Reinforced sand, Is a superposition of a sand model 
and an elastic "bilateral" model of the reinforcement. The "unilateral" TexSol 
model referred to §3.3 Is denoted Texsol. A particular model Is added denoted 
Spherical Texsol corresponding to the previous one with jiyj — 0. 

• Elasticity : Es = 200000 MPa, Us = 0.4, Ey, = 100000 MPa, = 0.3 

• Plasticity : Co = 50 MPa, Of ^0.1, H^^ 100 MPa, Hi = 100 MPa 
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Fig. 7. Zoom on the elastic range of the models 

The Sand and the Reinforced sand appear clearly as two elastic bounds for 
TexSol models (cf. figure 7). At this stage, the Spherical Texsol does not 
differ from the Sand. On the contrary, the Texsol is close to the Sand in 
compression and close to the Reinforced sand in traction. For the two loadings 
the limit models reveal to be the upper bounds. For a loading-unloading 
traction process, the Texsol model behaves almost like the Reinforced sand as 
expected (cf. figure 8). The Spherical Texsol does not improve significantly the 
Sand (cf. figure 8 and 9). Consequently, the Spherical Texsol does not account 
for the numerical results given in the figure 3 for the same kind of experiment 
- even roughly. 

4-3 Cyclic loading 

TexSol embankments may be submitted to vibrating solicitations. A cyclic 
test based on the test represented in the figure 10 (where the displacements 
are fixed on the lower side and the soUicitation managed by force) is performed 
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Fig. 8. Different material behaviours on a traction patch test 
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Fig. 9. Different material behaviours on a compression patch test 



to underline the contribution of the "unilateral" reinforcement due to the wire 
network. Some material parameters are changed to apply a greater amplitude 
of loading on it : Co = 150 MPa, Of = 0.02, Hk = 5 MPa, Hi = 1 MPa. 
Reinforced sand, Texsol and Sand are compared in the figures 11, 12 and 13. 
For the three models, the response tends to be stabilized after 20 cycles. 
But for Reinforced sand and Texsol the stabilization is reached before 10 
loops. Moreover the residual displacement of Texsol is 30 percents bigger than 
that of Reinforced sand and five times smaller than that of Sand. This last 
result highlights the advantages of TexSol reinforcement. An other effect of 
the "unilateral" wire in the Texsol model is clearly illustrated by the curvature 
changes when the displacement switches sign in the figure 12. 
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Fig. 10. Force managed traction / compression test 
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Fig. 11. Reinforced sand behaviour on a cyclic test (20 loops) 
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Fig. 12. Texsol behaviour on a cyclic test (20 loops) 
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Fig. 13. Sand behaviour on a cyclic test (20 loops) 
4-4 Compression test 

In soil mechanics it is usual to carry out a triaxial test with a prescribed 
confinement pressure (cf. fig 14). Considering the previous numerical results 
of the §4.2, the Spherical Texsol model is no more studied. Only the three 
other cases are compared in a loading compression test (the bulk mesh is 
described in [10]). The contribution of the wire in TexSol to the mechanical 




Fig. 14. Compression test 

strenght is illustrated by the spacial distribution of two stresses : the full stress 
(T and the wire stress cr^. The distribution of the full stress is identical in the 
three models with a level for Texsol between the two others. The main part of 
stress is located in the center of the bulk expect a localized concentration on 
the right lower corner. The contribution of the wire in the stress tensor (cr„^) 
is split into its deviatoric part and its spherical one (pressure). Both parts 
are identically null for Sand (cf. figure 16c and 17c). The elasticity of the 
reinforcement is activated only in tensile directions for the Texsol and in all 
directions for the Reinforced sand ; this explains the different full stress levels 
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a. Reinforced sand 



b. Texsol 



c. Sand 




Fig. 15. Equivalent Von-Mises stress 



a. Reinforced sand 



b. Texsol 
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Fig. 16. Wire equivalent Von-Mises stress 



in the figure 15 and the different wire stress levels in the figure 16. The nature 




Fig. 17. Wire pressure 



of the reinforcement due to the wire is clearly illustrated in the figure 17. The 

wire pressure in the Sand sample is of course identically null. It is negative in 
the Texsol wire (traction behaviour) according to the unilaterality condition 
expressed in the equation (20) whereas the pressure in the reinforcement of 
the Reinforced sand is almost everywhere positive. 
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5 Conclusion and prospects 

In this paper a coherent thermodynamical model is proposed to account for 

numerical experiments (because of the lack of real experiments on the TexSol). 
The keypoint is a "unilateral" elasticity which model the wire network. An 
elastic plastic model is superposed to the previous one to obtain both strain 
formulation and stress formulation when it is possible. Using a finite element 
method, we roughly validate the expected behaviour. 

The main perspective of this work is the identification of the mechanical 
parameters of the superposed model by series of numerical experiments in 
progress. In a more general framework orthotropic model is generally usefuU 
to model the wire network. 

The free energy ip considered in this work is postulated and in some cases we 
can write the free enthalpy via the Legendre-Fenchel transformation. Another 
approach should be to postulate the free enthalpy using a form similar to (19), 

The link between '^^y^ and ip'^ is an open question because in a three dimen- 
sional case the convexity of ip cannot be proved. 
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